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Abstract 

Populations of Castilleja austromontana (Scrophulariaceae) occur at high elevations (>2050 m) on the 
Madrean sky islands of southern Arizona, adjacent New Mexico and northern Mexico. These mountain- 
top habitat islands are cooler and wetter than the desert and semi-desert grasslands below. We examined 
41 vegetative, floral and seed characters for plants from six sky islands in Arizona. Univariate and 
multivariate analyses indicate strong morphological differentiation among populations. Discriminant func¬ 
tions analysis placed 88% of plants correctly, indicating that plants from different populations are mor¬ 
phologically distinct. Selection may contribute to the observed differentiation, but we have no basis for 
inferring the direction or strength of selection on vegetative traits. Sky island populations apparently share 
the same pollinators, so divergent selection on floral traits seems unlikely. Past and present gene flow 
may also be involved as our results conform, at least in part, to an isolation by distance model. The 
magnitude of character differences combined with inferences of population sizes suggests that genetic 
drift may also contribute to the differentiation observed. 


Evolutionary biologists have long recognized 
that islands, both real and habitat, provide special 
opportunities to study population differentiation 
and speciation due to geographic isolation (Mac- 
Arthur and Wilson 1967; Carlquist 1974; Grant 
1986, 1998 and papers therein; Wagner and Funk 
1995 and papers therein; Francisco-Ortega et al. 
1996; Givnish 1998). Gene flow will be reduced 
among populations on geographically isolated is¬ 
lands compared to contiguous populations. When 
gene flow is absent or rare, populations may dif¬ 
ferentiate over time as a result of selection or ge¬ 
netic drift (Wright 1943; Mayr 1954). Thus, island 
systems are ideal for studies of both the processes 
that hold populations together (i.e., gene flow) and 
those that result in differentiation (i.e., selection, 
genetic drift). 

In the Basin and Range province of the south¬ 
western United States and northwestern Mexico, a 
series of mountain ranges provide island-like hab¬ 
itats for plants and animals that require more mesic 
and temperate conditions than those offered by the 
deserts, scrublands and grasslands that now occupy 
the intervening basins (Marshall 1957). These Mad¬ 
rean sky islands developed as a series of deforma¬ 
tion mountains produced by continental rifting that 
occurred between 12 and 6 million years BP 
(Nations and Stump 1981; Morrison 1991). At pres¬ 
ent, 17 Madrean sky islands are high enough to 
support mixed conifer forest (ca. 2500 m elev.. 
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somewhat lower toward the north and higher to¬ 
ward the south). They form an archipelago that 
connects extensive and continuous highlands to the 
north (the Central Mountain Province and the Col¬ 
orado Plateau above the Mogollon Rim in Arizona) 
and south (the Sierra Madre Occidental that extends 
to within —75 km of the Arizona-Mexico border) 
(Fig. 1). As first documented by Shreve (1915; see 
also Heald 1950; Niering and Lowe 1984; Brown 
1994; Warshall 1995), the Madrean sky islands sup¬ 
port almost the same range of vegetation zones as 
occurs in all of the western United States, from 
open juniper woodlands just above semi-desert 
grasslands to spruce-fir forest on the summits of the 
highest mountains. 

Here we report morphological differentiation 
among populations of Castilleja austromontana 
Standi. & Blumer (Scrophulariaceae), an herba¬ 
ceous perennial that occurs on most of the Madrean 
sky islands. Castilleja austromontana is found 
above 2050 m elev., and therefore occurs as isolat¬ 
ed populations in the highest of sky island habitats. 
We used data from vegetative, floral and seed char¬ 
acters for individuals from six sky island popula¬ 
tions to address three questions: (1) Are populations 
distinguishable by the examined morphological 
traits? (2) If so, which traits most strongly differ¬ 
entiate the populations? (3) Can morphological di¬ 
vergence among populations be related to geogra- 
phy? 

Materials and Methods 

Castilleja austromontana is found as far north as 
the White Mountains of Arizona (Apache County) 
and the Sangre de Cristo Mountains of northern 
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Table 1. Location of populations of Castilleja austromontana from which samples were taken. Abbreviations 

following the name of each mountain range are those used in Tables 2 and 4. 

Catalina Mountains (Cat): Mount Lemmon, ca. 2 miles above Ski Valley, 2750 m elev., Pima County, Arizona. S. 
Slentz 44 (ARIZ). 

Santa Rita Mountains (SR): About 200 meters below Baldy Saddle, Old Baldy Trail, Mount Wrightson, 2500 m 
elev., Santa Cruz County, Arizona. S. Slentz and L. McDade 40 (ARIZ). 

Pinaleno Mountains (Pin): Hospital Flat, 2750 m elev., Graham Country, Arizona. S. McLaughlin and J. Bowers 
6431 (ARIZ). 

Rincon Mountains (Rin): Italian Ranch Trail, Italian Spring, Saguaro National Park East, 2400 m elev., Pima Coun¬ 
try, Arizona. S. Slentz 48 (ARIZ). 

Huaclmca Mountains (Hua): Carr Canyon Trail, 0.75 miles above Reef Area Campground, 2400 m elev., Cochise 
Country, Arizona. S. Slentz 43 (ARIZ). 

White Mountains (Wht): Lee Valley Recreation Area, ca. 3 miles south of Sunrise Lake on State Rt. 273, 2800 m 
elev., Apache County, Arizona. S. Slentz 45 (ARIZ). 


New Mexico and at least as far south as the limit 
of the sky island region in the northern Mexican 
states of Sonora and Chihuahua. Nothing is known 
about the age of the species nor is there a phylo¬ 
genetic hypothesis for Castilleja. Flowering of Cas¬ 
tilleja austromontana coincides with the summer 
wet season (July-September), and plants die back 
to the ground with the onset of winter. Our obser¬ 
vations (44 hours on the Huachuca, Catalina and 
White Mountains) indicate that flowers of C. aus¬ 
tromontana are visited and likely pollinated by 
hummingbirds. Individuals of Selasphorus rufus 
(Rufous Hummingbird) and possibly S. platycercus 
(Broadtailed Hummingbird) were observed to con¬ 
tact anthers and stigmas, and patches of pollen were 
seen on the heads of some visiting birds. There 
have been no detailed studies of the breeding sys¬ 
tems of plants of this species. The seeds of C. aus¬ 
tromontana have no special dispersal mechanism 
and seem to fall passively from dehiscing fruits. All 
Castilleja studied to date are autotrophic root par¬ 
asites (Heckard 1962, Marvier 1996); nothing is 
known about the nutritional status or hosts of C. 
austromontana. 

Collection Sites and Protocol. Field sites were 
located from specimens at the University of Ari¬ 
zona herbarium (ARIZ). We collected plants from 
one area on each of six different mountain ranges 
(Table 1; Fig. 1). During 1995 and 1996, we made 
at least one collecting trip to each mountain range 
during the flowering season to gather specimens for 
scoring vegetative and floral characters. These were 
followed by a later trip to collect mature fruits. In¬ 
dividual flowering plants were sampled by first de¬ 
termining the extent of the population, then walking 
a transect along the longest dimension of the pop¬ 
ulation and collecting from individuals at intervals 
chosen to achieve the desired sample size (15-20 
individuals per population) across this longest di¬ 
mension. Individual plants produce one to many 
stems from a root crown and are clearly distinct 
from other individuals. One flowering stem from 
each individual was snapped off just above the root 
stock at ground level. Only fresh flowers and their 
subtending bracts were scored (freshness was as¬ 


sessed by presence of pollen on recently dehisced 
anthers). 

Plants sampled for seed collection were chosen 
in the same manner. Infructescences from six to 11 
individuals were collected from each population ex¬ 
cept the White Mountains, where mature fruits 
could not be collected. Seeds were sampled only 
from mature, partially dehisced capsules; only 
filled, apparently viable seeds were used. Because 
individuals from which flowering stems were col¬ 
lected were not marked, it was not possible to col¬ 
lect infructescences from the same individuals. 

Morphometries. A total of 41 characters was re¬ 
corded from each flowering stem (35 measurements 
+ 6 meristics); of these, 11 were vegetative and 30 
were floral characters (Fig. 2, Appendix 1). Data 
were recorded from four leaves, three bracts, and 
three flowers per individual (Fig. 2). The sampling 
scheme for leaves and bracts specifically sought to 
capture positional differences in these structures: 
leaf and bract size and shape may vary along the 
stem or rachis, respectively, in Castilleja. We there¬ 
fore included the measurements for each individual 
leaf and bract in the analyses. We calculated means 
of corolla and calyx measurements for the three 
flowers sampled from each plant and used these 
means in all analyses. In addition, to examine leaf 
and bract shape, we combined length and width 
measurements of individual leaves and bracts. The 
ratio of corolla to galea length varies in Castilleja ; 
we determined this ratio for individual flowers and 
then calculated a mean value for each individual. 
The data set for univariate analysis thus includes a 
total of 39 characters (Table 2); of these, eight are 
ratios and the six calyx and corolla characters each 
represent means of three measurements. All mea¬ 
surements were log-transformed. Because entire 
bracts have cleft depths of zero, we added 0.5 to 
each value for depth of bract clefts prior to log- 
transformation. Meristic variables were square-root 
transformed. 

The intricate surface characteristics of the seeds, 
as viewed with the SEM, have proven useful in 
taxonomic work with Castilleja (Chuang and Heck¬ 
ard 1992; Lockwood 1992). Initial experiments 
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Fig. 1. The Madrean sky islands in southeastern Arizona and adjacent western New Mexico and northern Mexico. 
1500 m contour lines outline mountain ranges. Shaded areas are >2500 m elevation, corresponding to approximate 
range of habitat of Castilleja austromontana. Circles indicate relative placement of plant populations by multidimen¬ 
sional scaling (the MDS result was scaled so that its dimensions matched those of the physical map; C = Catalina 
Mountains, H = Huachuca Mountains, P = Pinaleno Mountains; R = Rincon Mountains; SR = Santa Rita Mountains; 
W = White Mountains). An X indicates the collection location for the White Mountain population. 


with techniques for preparing seeds (including de¬ 
hydration and critical point drying) indicated that 
there were no differences between seeds prepared 
using these methods versus simply placed directly 
onto SEM specimen stubs. Several seeds from the 
same fruit were then placed onto a stub; stubs were 
coated with 30-60 nm of gold-palladium and 
viewed using an l.S.l. DS-130 scanning electron 
microscope. Several photographs were taken of at 
least one seed from each of 6-11 plants per popu¬ 


lation for examination of seed characters. To facil¬ 
itate standardization of measurements, seeds chosen 
to be photographed were in the same orientation, 
i.e., with the longest dimension of the seed parallel 
to the long axis of the photo. The SEM was peri¬ 
odically recalibrated by photographing a microme¬ 
ter and adjusting the scale that appears on each 
photograph accordingly. Four measurements plus 
one meristic character were scored for each seed 
(Fig. 3, Appendix 1). To examine seed shape, we 
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Fig. 2. Castillejci austromontana. A. Schematic illustration of plant with enlarged terminal raceme; B. leaf; C. floral 
bract; D. calyx; E. corolla. Numbered dimensions indicate measurements taken (Appendix 1, Table 1): 1 Raceme 
length; 2 Leaf length; 3 Leaf width; 4 Bract length; 5 Bract width; 6 Anterior calyx cleft depth; 7 Lateral calyx cleft 
depth; 8 Posterior calyx cleft depth; 9 Galea length; 10 Corolla length. Drawn by S. Slentz. 


combined length and width measurements of indi¬ 
vidual seeds. When more than one seed was mea¬ 
sured from an individual plant, within-plant means 
of each character were used in the analyses. As 
with the vegetative + floral data, measurements 
were log-transformed and meristic data were 
square-root transformed. 

Statistical Analyses. The vegetative + floral data 
set was treated separately from the seed data be¬ 
cause the plants from which seeds were taken were 
not matched with those from which flowering stems 
were taken. All analyses were conducted in SAS 
(SAS Institute 1988) except that NTSYS (Rohlf 
1993) was used for multidimensional scaling (see 
below). One-way analyses of variance were used to 
determine which variables showed significant vari¬ 


ation among populations. Student-Newman-Keuls a 
posteriori tests were used to assess which popula¬ 
tions were significantly different from one another 
for each character. The vegetative + floral data set 
was subjected to principal components analysis 
(PCA) to examine relationships among characters. 
Ratios were omitted from the PCA because covari¬ 
ation between characters will be detected by mul¬ 
tivariate methods automatically. Because only five 
seed variables were measured, the seed data were 
not analyzed using PCA. Correlation matrices were 
used for PCA because both metric and meristic 
characters were included. 

Discriminant functions analysis (DFA) was used 
to determine how effectively the character data dis¬ 
tinguish plants from different mountain ranges. 
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Table 2. Results of ANOVA for Morphological Characters Ordered by Plant Structure. Results of the Stu- 
dent-Newman-Keuls a posteriori test are reported for each variable, with populations ordered from high values to low. 
See Table 1 for key to abbreviations of mountains. Lines below mountain abbreviations connect populations that are 
not significantly different. Asterisks indicate p-values that are judged significant by the sequential Bonferroni test. 


Character 


Results of ANOVA 

a posteriori test 

Vegetative : 








Leaf length 1 

F = 

1.83; 

5, 

83 

df; p 

= 0.1166 


Leaf length 2 

F = 

2.71; 

5, 

83 

df; p 

= 0.0255 

Wht Pin Hua Cat SR Rin 

Leaf length 3 

F = 

3.64; 

5, 

83 

df; p 

= 0.0050 

Wht Rin Pin Hua SR Cat 

Leaf length 4 

F = 

1.95; 

5, 

82 

df; p 

= 0.0950 


Leaf width 1 

F = 

8.51; 

5, 

83 

df; p 

< 0.0001* 

Pin Hua Rin SR Cat Wht 

Leaf width 2 

F = 

6.04; 

5, 

83 

df; p 

< 0.0001* 

Pin Hua Cat SR Rin Wht 

Leaf width 3 

F = 

2.01; 

5, 

83 

df; p 

= 0.0858 


Leaf width 4 

F = 

1.92; 

5, 

82 

df; p 

= 0.0999 


Leaf length/width 1 

F = 

0.48; 

5, 

83 

df; p 

= 0.7901 


Leaf length/width 2 

F = 

5.45; 

5, 

83 

df; p 

= 0.0002* 

Wht Cat Hua SR Rin Pin 

Leaf length/width 3 

F = 

6.12; 

5, 

83 

df; p 

< 0.0001* 

Wht Rin Hua Cat SR Pin 

Leaf length/width 4 

F = 

3.78; 

5, 

82 

df; p 

= 0.0039 

Wht Rin Hua Cat Pin SR 

No. lateral branches 

F = 

4.84; 

5, 

49 

df; p 

= 0.0011* 

Cat Rin Hua SR Pin Wht 

Apical stem width 

F = 

4.51; 

5, 

77 

df; p 

= 0.0012* 

Rin SR Hua Pin Wht Cat 

Basal stem width 

F = 

3.12; 

5, 

78 

df; p 

= 0.0129 

Hua Cat Rin Pin Wht SR 

Reproductive : 








No. of racemes 

F = 

3.64; 

5, 

63 

df; p 

• 0.0059 

Rin Cat Hua Pin SR Wht 

No. of flowers 

F = 

1.95; 

5, 

82 

df; p 

= 0.0943 


Raceme length 

F = 

1.65; 

5, 

82 

df; p 

= 0.1572 


Bracts: 








Bract length 1 

F = 

2.34; 

5, 

83 

df; p 

= 0.0489 

Wht SR Rin Hua Pin Cat 

Bract length 2 

F = 

2.74; 

5, 

83 

df; p 

= 0.0244 

SR Wht Pin Rin Hua Cat 

Bract length 3 

F = 

2.95; 

5, 

83 

df; p 

= 0.0168 

Wht SR Pin Rin Hua Cat 

Bract width 1 

F = 

1.20; 

5, 

83 

df; p 

= 0.3185 


Bract width 2 

F = 

1.62; 

5, 

83 

df; p 

= 0.1646 


Bract width 3 

F - 

2.56; 

5, 

83 

df; p 

= 0.0332 

SR Wht Hua Rin Cat Pin 

Length/Width 1 

F = 

2.57; 

5, 

83 

df; p 

= 0.0325 

Wht Rin Hua Pin SR Cat 

Length/Width 2 

F = 

2.09; 

5, 

83 

df; p 

= 0.0753 


Length/Width 3 

F = 

3.88; 

5, 

83 

df; p 

= 0.0033 

Pin Wht Rin Cat Hua SR 

Mean cleft depth 1 

F = 

4.49; 

5, 

83 

df; p 

= 0.0011* 

Pin Cat Wht Hua SR Rin 

Mean cleft depth 2 

F = 

4.31; 

5, 

83 

df; p 

= 0.0016 

Wht Pin Hua Cat SR Rin 

Mean cleft depth 3 

F = 

4.49; 

5, 

83 

df; p 

= 0.0011* 

Wht Pin Hua Cat SR Rin 

No. of clefts 1 

F = 

3.41; 

5, 

83 

df; p 

= 0.0075 

Pin Cat Wht SR Hua Rin 

No. of clefts 2 

F = 

2.26; 

5, 

83 

df; p 

= 0.0560 


No. of clefts 3 

F = 

3.36; 

5, 

83 

df; p 

= 0.0082 

Wht Pin SR Cat Hua Rin 

Calyx: 








Mean anterior cleft 

F = 

7.11; 

5, 

83 

df; p 

= 0.0001* 

Wht Pin Hua SR Rin Cat 

Mean lateral cleft 

F = 

6.29; 

5, 

83 

df; p 

= 0.0001* 

Pin Wht SR Hua Cat Rin 

Mean posterior cleft 

F = 

9.80; 

5, 

83 

df; p 

1 0.0001* 

Wht Pin SR Hua Rin Cat 

Corolla: 








Mean corolla length 

F = 

8.62; 

5, 

83 

df; p 

= 0.0001* 

Pin Hua Wht SR Rin Cat 

Mean galea length 

F = 

3.21; 

5, 

83 

df; p 

= 0.0106 

Pin Wht Hua Rin SR Cat 

Mean corolla/galea 

F = 

3.56; 

5, 

83 

df; p 

= 0.0057 

SR Pin Hua Wht Cat Rin 

Seeds : 








Seed length 

F = 

18.75 

i; 4, 37 df; p < 0.0001* 

Pin Hua Cat Rin SR 

Seed width 

F = 

4.52; 

4, 

37 df; p 

= 0.0022 

Pin Cat Hua Rin SR 

Seed length/width 

F = 

1.10; 4, 

37 df; p 

= 0.3724 


No. midline cells 

F - 

29.09 

4 

37 df; p < 0.0001* 

Pin Cat SR Hua Rin 

No. layers of testa 

F = 

15.28; 4 

37 df; p < 0.0001* 

Cat Hua Pin SR Rin 

Wall width 

F = 

12.13 

i; 4 

37 df; p < 0.0001* 

Pin Hua Rin Cat SR 
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Fig. 3. Scanning electron micrograph of seed of Ccistilleja austromontana ; seed length, width, and wall width were 
measured as indicated. Wall width was measured on a cell perpendicular to the line of vision. Number of midline cells 
= 14 on this seed as indicated by arrows. Enlarged cell in box illustrates layers of the radial walls (three in this case). 


Character data were log transformed only if their 
distributions were non-normal as DFA is not af¬ 
fected by scaling of individual variables (Manly 
1994). For seven plants from the Pinalenos, basal 
and apical stem widths were not measured; as a 
result, these plants were excluded by DFA. To max¬ 
imize the number of plants included, we conducted 
an additional DFA excluding the two characters that 
were missing for these plants. 

Non-metric multidimensional scaling (MDS) cre¬ 
ates map-like diagrams that optimize the placement 
of each population relative to all others; the algo¬ 
rithm minimizes “stress,” measured as the differ¬ 
ence between pairwise distances in the diagram 
produced by MDS and the original distances. MDS 
is widely used in the social and biological sciences; 
a number of authors have used this technique to 
study the relationship between genetic or phenotyp¬ 
ic distance and geographic distance (e.g., Wright 
and Ladiges 1997; Krauss 1996; Ruiz-Garcia 
1997). Guiller, Bellido, and Madec (1998) found 
MDS to be superior to other ordination methods; 


these authors provide a useful overview of the 
method and an explanation of its performance. Here 
we use MDS as a heuristic technique to visualize 
patterns of phenotypic distance in the context of 
geographic distances. Pairwise Mahalanobis dis¬ 
tances between populations were calculated in SAS 
and entered into NTSYS for the MDS analysis for 
each of the two Castilleja data sets. We scaled the 
diagram produced by this analysis so that its di¬ 
mensions matched those of the physical map de¬ 
picted in Fig. 1. The MDS results were then over¬ 
laid on the geographic map of the sky island region 
to compare relative phenotypic and geographic dis¬ 
tances (Fig. 1). 

Results 

Twenty-eight of 39 vegetative + floral characters 
and all seed characters except ratio of seed length/ 
width varied significantly among populations (Ta¬ 
ble 2). Twelve vegetative + floral characters and 
four seed characters vary significantly among pop- 
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Table 3. Eigenvectors on First Three Principal Components for Characters Included in Principal Component 
Analysis of the Vegetative and Reproductive Data Excluding Seeds and Ratios. Analysis included all characters. 
The first three principal components account for 52.9% of the total variation. See Appendix 1 for explanation of 
characters. 


Characters 

PCI 

PC2 

PC3 

Bract length 2 

0.2736 

0.1076 

-0.1530 

Leaf length 4 

0.2703 

0.0622 

-0.0533 

Mean anterior calyx cleft depth 

0.2653 

-0.0975 

-0.0242 

Leaf length 3 

0.2630 

0.0495 

-0.0786 

Bract length 1 

0.2619 

0.1452 

-0.1370 

Mean posterior calyx cleft depth 

0.2528 

-0.1973 

0.0011 

Mean corolla length 

0.2461 

0.0168 

-0.0076 

Bract length 3 

0.2402 

-0.0080 

-0.2036 

Bract width 2 

0.2344 

-0.0382 

0.1694 

Leaf width 3 

0.2262 

0.0857 

0.1295 

Leaf width 4 

0.2188 

0.1475 

0.1563 

Mean galea length 

0.1916 

0.0443 

-0.0159 

Raceme length 

0.1907 

0.0358 

0.1193 

Number of racemes 

-0.1825 

0.1042 

0.2325 

Number of lateral branches 

-0.1807 

0.1228 

0.2234 

Bract width 1 

0.1739 

0.2196 

0.1607 

Bract width 3 

0.1731 

0.1393 

0.0733 

Mean lateral calyx cleft depth 

0.1609 

-0.1291 

-0.0052 

Leaf length 2 

0.1422 

0.0328 

0.0363 

Number of flowers per raceme 

0.1244 

-0.0002 

0.2762 

Bract 3, mean cleft depth 

0.1225 

-0.2549 

0.1069 

Bract 2, mean cleft depth 

0.1204 

-0.3384 

0.1773 

Leaf width 2 

0.1194 

0.2331 

0.1976 

Bract 2, number of clefts 

0.1027 

-0.2743 

0.2194 

Bract 3, number of clefts 

0.0592 

-0.3596 

0.2684 

Leaf length 1 

-0.0502 

0.2115 

0.2808 

Basal stem width 

-0.0415 

0.1710 

0.3017 

Leaf width 1 

-0.0305 

0.3290 

0.2622 

Bract 1, number of clefts 

-0.0295 

-0.2443 

0.2698 

Bract 1, mean cleft depth 

0.0266 

-0.2354 

0.2871 

Apical stem width 

-0.0262 

0.1100 

0.1219 

Eigenvalue 

8.95 

4.14 

3.32 

% variance explained 

28.88% 

13.36% 

10.69% 


illations even after adjusting threshold probability 
values via the sequential Bonferroni test (Rice 
1989) (Table 2). Characters that vary significantly 
among populations are from diverse plant structures 
(Table 2). Leaf shape (i.e., length/width) and size 
differ among populations, with plants from the 
White Mountains having leaves that are relatively 
longer and narrower than those of plants from the 
other sky islands. Plants from the White Mountains 
also have relatively few branches. Bracts differ in 
cleftedness (i.e., both number and depth of clefts), 
as well as in shape (length/width) and length among 
sky islands. Plants from the White Mountains and 
the Pinalenos have calyces that are more deeply 
cleft than those of other sky island populations. Co¬ 
rollas differ in total length and in galea length; 
flowers of plants from the Pinalenos are largest 
whereas those from the Catalinas are smallest. 
Plants from the Pinalenos have long, wide seeds 
with many cells along the midline and wide walls 
compared to seeds from plants from the Santa Ritas 
or Rincons. Seeds from plants from the Catalinas 
have the most layers of the radial wall and those 


from the Rincons the fewest; seeds from plants in 
the Pinalenos are intermediate for this character. 

For the analysis including vegetative + floral 
data, the first three principal components each ac¬ 
count for >10% of the variation among individuals, 
but together explain only 52.9% of total variation. 
PCI accounts for 28.9% of the standardized vari¬ 
ance in the data. The characters that contribute most 
to PCI are length of bract 2, depths of the calyx 
clefts, and corolla length (Table 3). Raceme length, 
bract length and width, and leaf length and width 
also contribute substantially. The characters con¬ 
tributing most to PC2 involve cleftedness of bracts 
(i.e., both number and depth of clefts), and bract 
length (Table 3). This component explains an ad¬ 
ditional 13.4% of the standardized variance. Char¬ 
acters contributing most to PC3 are basal and apical 
stem width, number of flowers per raceme, and leaf 
width. The third principal component accounts for 
10.7% of the standardized variance. 

The next two principal components each account 
for more than 5% of the total variation for a cu¬ 
mulative total of 68% of total variation accounted 
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Table 4. Number (Above) and Proportion (Below) of 
Plants From Source Locations (Left) That Were Iden¬ 
tified by Discriminant Functions Analysis of Vegeta¬ 
tive + Floral Data as From the Locations Listed 
Above. Characters were excluded for which data were 
missing for some plants. See Table 1 for key to abbrevi¬ 
ations of mountains. 



Cat 

Hua 

Pin 

Rin 

SR 

Wht 

Cat 

10 

1 

0 

2 

0 

0 


76.9 

7.7 

0.0 

15.4 

0.0 

0.0 

Hua 

1 

12 

0 

0 

0 

0 


7.7 

92.3 

0.0 

0.0 

0.0 

0.0 

Pin 

0 

0 

15 

0 

0 

0 


0.0 

0.0 

100.0 

0.0 

0.0 

0.0 

Rin 

1 

0 

0 

8 

2 

0 


9.1 

0.0 

0.0 

72.7 

18.2 

0.0 

SR 

1 

0 

0 

1 

16 

0 


5.6 

0.0 

0.0 

5.6 

88.9 

0.0 

Wht 

0 

0 

1 

0 

0 

15 


0.0 

0.0 

6.3 

0.0 

0.0 

93.8 


for by the first five principal components. Sixteen 
principal components must be considered in order 
to explain more than 95% of the variation. Follow¬ 
ing Legendre and Legendre (1998) and considering 
only PCs with eigenvalues >1 restricts us to the 
first eight PCs which account for only 81% of the 
variation. Although there are no established guide¬ 
lines for judging the power of a PC A to explain 
variation, Sneath and Sokal (1973) suggest that 
three PCs satisfactorily illustrate patterns of varia¬ 
tion in most systematic studies. This is not the case 
here, indicating that many characters vary indepen¬ 
dently of others such that morphological variation 
among sky island populations of Castilleja defies 
reduction to a few orthogonal principal compo¬ 
nents. 

Discriminant functions analysis including all 
plants (i.e., omitting basal and apical stem width 
which were missing for seven plants from the Pin¬ 
alenos) indicated that plants can be linked to their 
source population with considerable accuracy: 88% 
(76 of 86 plants) are correctly assigned (Table 4). 
Nine of ten errors involved plants from the southern 
four populations (i.e., the Catalinas, Rincons, Santa 
Ritas and Huachucas) and no plants from the Pin¬ 
alenos were misplaced. 

Seeds of C. austromontana are of Lockwood’s 
(1992) Type 1; the outer membrane is essentially 
completely ruptured and lost at maturity, the outer 
tangential wall is smooth, and the radial walls of 
the “cells” are unornamented. Seed characters vary 
among populations (Table 2) and discriminant func¬ 
tions placed more than 3 A of plants correctly based 
on seed characters (Table 5). Errors were most fre¬ 
quent regarding seeds from the Huachucas; nearly 
half were incorrectly placed. However, no seeds 
from the Pinalenos were misclassified as being 
from other sky islands nor were seeds from other 
populations placed on the Pinalenos. Seed mor- 


Table 5. Number (Above) and Proportion (Below) of 
Plants From Source Locations (Left) That Were Iden¬ 
tified by Discriminant Functions Analysis of Seed 
Data as From the Locations Listed Above. See Table 
1 for key to abbreviations of mountains. 



CAT 

HUA 

PIN 

RIN 

SR 

CAT 

5 

1 

0 

1 

1 


62.5 

12.5 

0.0 

12.5 

12.5 

HUA 

1 

6 

0 

2 

2 


9.1 

54.6 

0.0 

18.2 

18.2 

PIN 

0 

0 

9 

0 

0 


0.0 

0.0 

100.0 

0.0 

0.0 

RIN 

0 

0 

0 

6 

0 


0.0 

0.0 

0.0 

100.0 

0.0 

SR 

2 

0 

0 

0 

6 


25.0 

0.0 

0.0 

0.0 

75.0 


phology among plants on the Pinalenos is distinc¬ 
tive; these seeds are large with wider walls than 
those of plants from other sky islands. 

Superimposing the diagram from the MDS anal¬ 
ysis of vegetative + floral characters on the map of 
the sky island region indicates that, among the Ca¬ 
talinas, Rincons, Santa Ritas, and Huachucas, phys¬ 
ical and relative phenotypic distances are roughly 
congruent (Fig. 1). Relative to these four, the White 
Mountains are placed nearer and Pinalenos farther 
by MDS based on phenotypic distances than their 
geographic distances suggest. 

The MDS analysis of seed data (not shown) 
placed the Pinaleno population distantly from the 
other four populations which were proximate and 
approximately equidistant from each other. 

Discussion 

Our results demonstrate morphological differen¬ 
tiation among populations of Castilleja austromon¬ 
tana on Madrean sky islands. DFA associated 88% 
of plants with the population from which they were 
collected, and PCA indicates that many vegetative 
and reproductive characters contribute to this dif¬ 
ferentiation. At one extreme, plants on the White 
and Pinaleno Mountains tend to be relatively small 
in stature, with long leaves and bracts, long corol¬ 
las, and deep bract and calyx clefts. At the other 
extreme, on the Catalina and Rincon Mountains, 
plants tend to be relatively tall with smaller leaves 
and flowers, and shallow clefts in bracts and caly¬ 
ces. Plants on the Huachuca and Santa Rita Moun¬ 
tains have structures that are intermediate in size 
and shape between these extremes. 

It is important to note that we do not know the 
cause of this morphological variation. It is possible 
that the characters are phenotypically plastic and 
that our results reflect environmental differences 
among mountains. These plants are likely root par¬ 
asites, so that common garden or reciprocal sowing 
or transplant approaches to examine this possibility 
will be difficult at best. Because the sizes and 
shapes of reproductive structures are generally un- 
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derstood to be under selection associated with 
achieving successful pollination and thus reproduc¬ 
tion, these are often less phenotypically plastic than 
vegetative characters (e.g., Bradshaw 1965; Davis 
1983; Apelgren and Lernstal 1991). Our results in¬ 
dicate that both vegetative and floral characters 
vary among sky island populations, suggesting that 
we are not dealing solely with phenotypic variants. 
However, data from other sources will be necessary 
to confirm that the differentiation that we see is 
genetically based. To the extent that the characters 
studied here have a genetic basis, they may also be 
genetically correlated and thus not reflect a suite of 
independent genetic differences. However, PC A 
suggests that many characters are independent of 
each other, at least in their phenotypic expression. 

Selection, gene flow, and genetic drift may all 
contribute to the patterns of differentiation that we 
have documented. Contrasting selective regimes 
may be responsible for the morphological differ¬ 
ences among plants from different sky island pop¬ 
ulations. We have no basis at present for specific 
hypotheses regarding the direction or magnitude of 
selection except that disruptive selection on corolla 
morphology is unlikely because pollinator relation¬ 
ships appear to be uniform. 

Some inferences regarding gene flow are possi¬ 
ble based upon geographic distances, pollinator re¬ 
lationships, and vegetation history. If gene flow 
(whether by pollen or migrants) is negatively relat¬ 
ed to geographic distance, then proximate popula¬ 
tions should be more similar than populations that 
are distant and thus rarely exchange genes. The 
physical distances among the Catalina, Rincon, 
Santa Rita, and Huachuca populations are fairly 
congruent with phenotypic distances among plant 
populations from these ranges (Fig. 1) in accord 
with a gene flow or isolation by distance explana¬ 
tion. In contrast, the placement by MDS based on 
phenotypic distance of the White Mountain and 
Pinaleno populations relative to the other four is not 
congruent with geographic distances (Fig. 1). In 
this light, the behavior of pollinators may be infor¬ 
mative. 

As described above, observations suggest that 
the same hummingbird species visit flowers of C. 
austromontana across the species’ range. Individ¬ 
uals of Selasphorus platycercus breed on the sky 
islands at the elevation of the Castilleja populations 
and thus should move pollen within ranges. In con¬ 
trast, Rufous-tailed Hummingbirds (S. rufiis) pass 
through these mountains during the late summer 
and early fall migration south to their winter range 
(the birds fly at lower elevation during their spring 
migration north because they precede the spring 
thaw at high elevation) (Grant and Grant 1967; 
Johnsgard 1983). If relative phenotypic distance re¬ 
flects gene flow, then our data suggest a more direct 
route between the southern populations and the 
White Mountains than would be expected based on 
geographic distance. Unfortunately, migration 


routes are understood only in broad geographic 
terms, not at the level of routes followed by indi¬ 
vidual birds. 

Vegetation history of these mountain ranges 
yields some ideas about past patterns of gene flow. 
During the late Pleistocene and Holocene, cycles of 
climate change were associated with glaciations 
farther north and at higher elevations. Fossil pack- 
rat middens provide direct evidence of vegetation 
changes associated with climatic changes during 
the last 40,000+ years (Van Devender, Thompson 
and Betancourt 1989; Betancourt et al. 1990 and 
references therein). During periods of glacial max¬ 
ima, vegetation zones on the sky islands were sub¬ 
stantially lower than they are today. The most com¬ 
plete evidence, from the most recent glacial maxi¬ 
mum (i.e., 20,000 to 8,000 years BP), indicates that 
the habitat in which C. austromontana occurs 
ranged 500-600 m lower than at present. This 
would presumably have increased the area of suit¬ 
able habitats on the sky islands and also increased 
connectivity among them. Higher connectivity 
would have likely increased gene flow between 
populations that are proximate (e.g., Catalinas and 
Rincons) and between those that are separated by 
valleys of relatively high elevation (e.g., Huachucas 
and Santa Ritas). Both the geographic distance and 
topographic complexity between these four sky is¬ 
lands and the Pinalenos and White Mountains pre¬ 
clude simple predictions of the impact of climate 
change on gene flow. 

Finally, genetic drift occurs in populations that 
are isolated, especially if they are small (Lande 
1976; Hard and Clark 1989). We do not have cen¬ 
sus data for these populations, but some inferences 
are possible based upon areal extent of suitable 
habitat. The area inhabited by C. austromontana in 
the White Mountains is vast such that this popula¬ 
tion should be least susceptible to chance effects. 
In contrast, the area of suitable habitat for C. aus¬ 
tromontana on all of the sky islands is far more 
limited and all of these populations may be subject 
to drift. This is especially true if warm interglacial 
periods produced an upward elevational shift in cli¬ 
matic zones comparable to the well-documented 
downward shift associated with glacial maxima. 

Lande (1976:32) derived an equation to estimate 
the population size (N*) below which drift may in¬ 
fluence morphological evolution: 

_ (1.96) 2 h 2 t 

(z/5) 2 

where h 2 is heritability, t is number of elapsed gen¬ 
erations, z is the difference in character means be¬ 
tween two samples (i.e., the magnitude of morpho¬ 
logical change observed), and 8 is the standard 
deviation of the character. We estimated values for 
this population size threshold in C. austromontana 
using our data for corolla length from the Huachu¬ 
cas and Santa Ritas. Heritability has not been esti¬ 
mated for corolla length in C. austromontana , but 
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Campbell (1996) obtained h 2 values between 0.24 
and 0.74 for another Asteridae, Ipomopsis aggre¬ 
gate! (Polemoniaceae). It seems likely that h 2 for 
corolla length in C. austromontana falls within this 
rather large range. Further, we can estimate time 
since these populations were in contact: because of 
the relatively high basin between these two ranges, 
we hypothesize that there would have been nearly 
continuous habitat for C. austromontana during the 
last glacial maximum, ca. 10,000 years BP. Gen¬ 
eration time is not known for these plants and we 
used 5 and 10 yr values as reasonable for herba¬ 
ceous perennials. The range of parameters explored 
gives estimates from 21,000 (h 2 = 0.24, generation 
time = 10 yr) to 130,000 individuals (h 2 = 0.74, 
generation time = 5 yr); that is, drift cannot be 
ruled out if population sizes are lower than these 
values. Based on field experience and on areas of 
suitable habitat among sky islands, it seems likely 
that only the White Mountains could have as many 
as 20,000 individuals. These results are consistent 
with the hypothesis that drift is at least partly re¬ 
sponsible for the differences that we see among sky 
island populations. 

There are no other quantitative studies of plant 
population differentiation among Madrean sky is¬ 
land plants. Maddison and McMahon (in press) 
have documented substantial morphological and 
behavioral differences among male spiders ( Ha- 
bronattus : Salticidae) on different sky islands. 
Snails of the genus Sonorella (Helminthoglyptidae; 
Miller 1967) and beetles (, Scaphinotus , Carabidae; 
Ball 1966) are reported to be highly differentiated 
among the Madrean sky islands although quantita¬ 
tive or phylogeographic studies are lacking. Such 
studies are underway for spiders (S. Masta personal 
communication) and lizards and snakes (K. Zamu- 
dio and H. Greene, Cornell University, personal 
communication). Other sky island systems (e.g., the 
tepuis of Venezuela) have been the focus of floristic 
and evolutionary studies above the species level 
(Givnish et al. 1997), but we are not aware of any 
research to document intraspecific patterns of dif¬ 
ferentiation there. Plant populations on true island 
systems have been studied by a number of research¬ 
ers, but patterns of intraspecific differentiation have 
received relatively little attention (but see Francis- 
co-Ortega et al. 1993; Hotta et al. 1985), in part 
because there appears to be little genetic divergence 
within many island lineages (Givnish 1998). 

In sum, our results confirm that the Madrean sky 
islands are likely to be evolutionary hot spots and 
that further research is warranted. In the case of C. 
austromontana , additional populations should be 
studied, particularly from the southern portion of 
the Madrean sky island province (i.e., the Mexican 
sky islands between the international border and the 
“mainland” of the Sierra Madre Occidental). Ge¬ 
netic data that can be ordered (e.g., DNA sequenc¬ 
es) would permit us to test more explicitly our ideas 
regarding gene flow and genetic drift, as well as to 
add a phylogeographic component (sensu Avise 


1989, 1994; Templeton 1998). Further, parallel 
studies on other groups of plants with both similar 
and contrasting life histories, pollination systems 
and dispersal mechanisms will deepen our under¬ 
standing of both the correlates and determinants of 
population differentiation. 
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Appendix 1. Explanation of characters measured or counted for Castilleja a ustromontana. See Figures 2 and 
3 for illustrations of structures measured. Terminology follows Chuang and Heckard (1992). 


Vegetative: 

Basal Stem Width was measured just above the root stock. 

Apical Stem Width was measured directly below the inflorescence. 

Leaf Length (Fig. 2B-2) and Leaf Width (Fig. 2B-3) were measured for four leaves per sampled stem including the 
most proximal (leaf 1), most distal (i.e., last vegetative leaf below the most proximal floral bract; leaf 4), and two 
intermediates at ca. Vi and % of total stem height. 

Number of Lateral Branches includes both vegetative and reproductive branches. 

Floral: 

Number of Racemes is number of branches bearing an inflorescence. 

Number of Flowers on the raceme selected for measurement of floral characters. 

Raceme Length was measured from directly below the most proximal floral bract to the tip of apical bud (Fig. 

2A-1). 

Bract, Calyx and Corolla: three nodes were sampled and scored for bract, calyx and corolla characters; these were 
the most proximal (1), most distal (3), and approximately median (2) flowering nodes. 

Bract Length (Fig. 2C-4) and Bract Width (Fig. 2C-5) 

Number of Clefts per Bract was counted on each of the three sampled bracts; bracts have from zero (entire) to four 
deep clefts on each lateral margin. 

Mean Cleft Depth per Bract of all clefts on a given bract. Because number and position of clefts varied among 
bracts, it was not possible to compare individual clefts. This variable attempts to capture ‘cleftedness,’ which 
clearly varied among plants, without comparing individual clefts. 

Calyx : Depths of the anterior (Fig. 2D-6), lateral (Fig. 2D-7), and posterior (Fig. 2D-8) calyx clefts were measured 
as the straight line distances from the base of each cleft to the apex of an adjacent calyx lobe. For analysis, 
measurements for the three sampled calyces were combined as mean depths of anterior, lateral and posterior clefts. 

Corolla Length is from the apex of the posterior lobes of the corolla to the posterior base of the corolla (Fig. 2E- 
10). For analysis, measurements for the three sampled corollas were combined as mean corolla length. 

Galea Length is from the apex of the posterior lobes of the corolla across the mouth of the corolla to the apex of 
the anterior lip (Fig. 2E-9). For analysis, measurements were combined as mean galea length. 

Seeds: All variables were scored from SEM photomicrographs (Fig. 3). When more than one seed per individual was 
photographed, measurements or counts were averaged for analysis. 

Seed Length was measured from the hilium to the opposite end of the long axis of the seed. 

Seed Width was measured at the widest point of the seed as photographed in the standard orientation. 

Number of Layers of the Radial Wall is the number of layers of “scalariform thickenings on the radial walls” 
(Lockwood, 1992:225). 

Number of Midline Cells : We counted number of cells in the longest row along the midline of the seed. 

Radial Wall Width: Width of the radial wall of testa cells measured at the widest point on a cell directly perpendicu¬ 
lar to the line of sight. 

Derived Variables: These ratios were calculated for each leaf, bract, and corolla measured; ratio of seed length/seed 
width was calculated based on within-plant mean length and width values (see Table 2). 

Leaf Length/Leaf Width: This variable describes leaf shape; leaves with low ratios are more ovate than those with 
higher ratios. 

Bract Length/Bract Width: This variable describes bract shape; bracts with low ratios are more ovate than those with 
higher ratios. 

Corolla length/Galea length: This variable describes the corolla shape. Flowers with low ratios have long galeas 
relative to corolla length, and those with higher ratios have short galeas relative to corolla length. 

Seed Length/Seed Width: This variable describes seed shape; seeds with low ratios are more ovate than those with 
higher ratios. 




